Signatures of Radiation Reaction in Ultra-Intense Laser Fields 
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We discuss radiation reaction effects on charges propagating in ultra-intense laser fields. Our analysis 
is based on an analytic solution of the Landau-Lifshitz equation. We suggest to measure radiation 
reaction in terms of a symmetry breaking parameter associated with the violation of null translation 
invariance in the direction opposite to the laser beam. As the Landau-Lifshitz equation is nonlinear 
the energy transfer within the pulse is rather sensitive to initial conditions. This is elucidated by 
comparing colliding and fixed target modes in electron laser collisions. 



I. INTRODUCTION 

The problem of classical radiation reaction (RR) has 
vexed generations ofphysicists since its first formulation 
in 1892 by Lorentz Jl|, |2| • Following important contribu- 
tions by Abraham [3| and others, the equation describing 
the back reaction of the radiation field on the motion of 
the radiating charge has been cast in its final covariant 
form by Dirac in 1938 Q- It is now aptly called the 
Lorentz- Abraham-Dirac (LAD) equation. The relevant 
body of literature has become enormous and we refer to 
the recent monographs 0, and, in particular, to the 
preprint (3] for an overview of the historical development 
and extensive lists of references. 

A particularly compact way of writing the LAD equa- 
tion, say for an electron (mass m, charge e) is 



mu = + To P mu , 



(1) 



where u denotes the electron 4- velocity, F = eFw/c the 
Lorentz 4-force in terms of the field strength tensor, F, of 
the externally prescribed field and dots derivatives with 
respect to proper time, r. The second term on the right 
is the RR force, i^RR, which is characterised by the ap- 
pearance of the time parameter 



TO ^ h-eic ~ 2 fm/c ~ 10 



-23 



(2) 



This is the time it takes light to traverse the classical 
electron radiufl 



e^/47rmc^ = aXc — 3fm , 



(3) 



or, from the second term, the electron Compton wave 
length reduced by a factor a ~ 1/137, the fine structure 
constant. Obviously, the time and length scales involved 
are typical for higher order QED corrections (or even 
strong interactions, i.e. QCD) - a hint that the classi- 
cal LAD equation ([1]) may not capture all the physics at 
these microscopic (and essentially quantum) scales. Fi- 
nally, the projection P = 1 — ui»u/<? in ([1]) guarantees 



that 4-acceleration and velocity are Minkowski orthogo- 
nal. This follows upon differentiating the on-shell condi- 
tion, vP' = c^, which, of course, is Einstein's postulate on 
the universality of the speed of light, c. As the velocity 
u is time-like, the acceleration ii is space-like, < 0. 



II. ESTIMATING RADIATION REACTION 

The LAD equation (jT]) is of third order in time deriva- 
tives (ii = 'x) and hence suffers from a number of patholo- 
gies such as runaway solutions and pre-acceleration. One 
way to overcome this is by iteration, assuming that 
-Frr, F - which amounts to working to first order in To. 
This in turn implies a 'reduction of order' in derivatives 
and results in the Landau-Lifshitz (LL) equation 



F + tqTF . 



(4) 



Hence, one replaces the offending 'jerk' term, mu, in 
fl]) by the proper time derivative of the Lorentz force 
[6| where the ii term is evaluated to lowest order in To, 
giving 



F ^ - Fu H W^u + 0(to) 

c mc^ 



(5) 



For alternative derivations of the LL equation resolving 
mathematical intricacies related to regularisations of the 
point particle concept we refer to O,!!!- 

The LL equation Q was derived under the assumption 
of a small reaction force, Frr <C F. Let us elucidate the 
physics involved somewhat further by assuming that the 
external field is produced by a laser described by a plane 
wave with light-like wave vector fc, fc^ — 0. An electron 
'approaching' the laser field with initial 4- velocity uo will, 
in its rest frame, 'see' a wave frequency given by the 
scalar product. 



k ■ uq . 



(6) 
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^ We employ Heaviside-Lorentz units with fine structure constant ^ We denote the tensor product U/^Ui, in index-free notation with 
a = e^/^Trhc = 1/137. the standard symbol, 
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At this point one has to distinguish between two points of 
view. If both k and uq are simuhaneously Lorentz trans- 
formed the frequency D,q remains, of course, invariant. 
On the other hand, one may think of fc, the wave vector 
of the laser as measured in the lab, as a distinguished 4- 
vector that breaks explicit Lorentz invariance (selecting 
a specific photon energy and beam direction). Different 
choices of initial conditions (i.e. uq) then characterise the 
relation between the initial rest frame of the electron and 
the lab frame. In what follows we will adopt this second 
point of view. 

The temporal gradients in ((S]) will be of the order of the 
laser period, F ~ Ogi^, so that the relative magnitude of 
the reaction force becomes 



Ft 



RR 



F 



(7) 



with the inequality required for the validity of the LL 
approximation. Consider now a head-on collision in the 
lab where the laser wave vector and electron velocity are 
given by 



fc = cj/c(l,z). Mo = 7oc(l, -/3oz) , 



(8) 



with the usual relativistic gamma factor, 70 = Ee/mc^ 
measuring the electron energy E^ in units of mc^. Such 
an electron then 'sees' a laser frequency that is Dopplcr 
upshifted according to 



^0 = 7o(l + = e''°uj ~ 2-fQUj 



(9) 



with the last identity holding for 70 3> 1. This boost in 
laser frequency is just the usual energy gain due to collid- 
ing versus fixed target mode (which, of course, are related 
by a longitudinal Lorentz boost with rapidity ^o)- If we 
define dimensionless photon energies in the co-moving 
and lab frames. 



1^0 



2 ' TT-i ^2 



the RR parameter r from ([7]) becomes 



r = \avQ ~ ^a-iQV ~ 10 ^701^ 



(10) 



(11) 



and hence is of order tq . The LL equation ([4]) would thus 
give us the radiated power to order Tq . If we follow our 
philosophy of neglecting second-order terms it suffices to 
express the acceleration via just the Lorentz force. 



F = eFu/c = eE 



(13) 



where wc have introduced the space-like 4- vector E cor- 
responding to the electric field 'seen' by the electron. If 
we now assume trans versality of our laser fields, k-E = 0, 
we immediately derive a very useful conservation law by 
dotting k into (fT3|) . 



Q. = k-u = Q 



(14) 



In other words, without RR, the electron always 'sees' 
the same laser frequency on its passage through the laser 
beam. 



Q, = k ■ u — k ■ Uq ^ 



(15) 



If we now plug the Lorentz equation (jl3)) into ((12)) the 
average energy loss per (reduced) laser period 1/fio in 



units of 1 



R 



becomes 



e^{E^) 



(16) 



Interestingly, we recover our RR parameter r from ^ 
and (jlip multiplying a new quantity, the dimensionless 
laser amplitude oq. This measures the energy gain of an 
electron traversing a laser wavelength, Xl = c/fJo, in a 
field of average strength (— iJ^)^/^ in units of mc^ . Obvi- 
ously, when this (purely classical) parameter becomes of 
order unity, the electron motion is relativistic. Note that 
ao is Lorentz invariant if k and Uq are transformed simul- 
taneously. Gauge invariance is shown [l3j by expressing 
it in terms of the field strength, F which is conveniently 
rendered dimensionless by introducing 



so that Oq finally becomes 



(mo, (F2)"o) 



(17) 



(18) 



For an optical laser v ~ 10^^ so that r ~ 10^^7o. Thus, 
to boost this to order unity (such that reaction equals 
Lorentz force) requires 70 — 10^, i.e. electron energies of 
order 10^ TeV. These can only be produced in gamma-ray 
bursts, but not (currently) in labs. The ground breaking 
laser pair production ("matter from light") experiment 
SLAG E-144, for instance, was utilising the 50 GeV SLAG 
linear collider implying 70 ~ 10^ and r ~ 10"'^ [l^ . 

The standard way of quantifying radiation by accel- 
erated charges is via Larmor's formula for the radiated 
power, the relativistic incarnation of which may be writ- 
ten as 



-niTQU > 



(12) 



The eiiergy loss parameter (|16p was previously employed 
in pH - [l7| . It suggests that for substantial radiation the 
small RR parameter r needs to be compensated by large 



values of 



The magnitude of oq is most easily es- 



timated by introducing Sauter's critical electrical field, 
Es = m^c^/eh = 1.3 x 10^^ V/m [H] and the associated 
intensity. Is = cEg = 4.2 x 10^^ W/cm^. For a given lab 
intensity / we then have 



ao 



y Is 



(19) 



In view of the current record intensity of / = 10^^ W/cm^ 
one can envisage oq values of about 10'^ for the not 
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too distant future [20|, |2l| . We have seen already in (fTTI) 
that large gamma factors (colliding mode) yield a fur- 
ther increase of radiative losses. In addition, the losses 
accumulate over successive laser periods. After, say, N 
cycles one expects a total relative change of the electron 
gamma factor given by 



= 2Tr NR = —Na-foiyal 
70 3 



(20) 



Thus, the smallness of av ~ 10~* may be compensated 
by pulse duration, N, initial electron energy, 70, and in- 
tensity, Qq. 



III. MODELLING THE LASER 

The simplest model of a laser (beam) is provided by 
a plane wave with a field tensor depending solely on the 
phase, F = F(0), cj) = k ■ x, and obeying transversality, 
Ffc = 0. If we choose fc as in (151 we have 



n 



k ■ u 
k ■ X 



oj{u° -u^)/c = uju-/c , (21) 
uj{x° - x^)/c = LUX- /c . (22) 



In other words, the laser field F only depends on the 
light- front or null coordinate, x^ = ct — z [22| . 

Plane waves are invariantly characterised as null fields 
[23I [2^ for which both scalar and pscudoscalar invariants 
vanish, 



itrF^ 
itrFF 



(23) 
(24) 



The vanishing of J?' implies that the energy momentum 
tensor of a plane wave is just F^, 



cT = ¥^ -yi = F^ 



(25) 



This is the only nontrivial power of field strength as F 
is nilpotent with index 3, i.e. F'^ = 0, which will be im- 
portant when we solve the equations of motion in such a 
field. 

We emphasize that there is no intrinsic invariant scale 
associated with a null field. The only way to associate 
a nonvanishing invariant is by using a probe (dubbed 
"third agent" in [lH) such as an electron or a (non-laser) 
photon. This naturally leads to the invariant amplitude 
ao as defined in psp which explicitly depends on the 
probe electron 4-velocity. Defining the energy density 
of the laser 'seen' by the electron as 



(26) 



we see that ([T8|) precisely represents the dimensionlcss 
version of this energy density. 



(27) 



Typically, the plane wave modelling the laser will be 
pulsed, i.e. of finite duration in phase (j). We accommo- 
date this situation by paramcterising the dimensionlcss 



field strength F as follows. We assume the plane wave 
field to be linearly polarised along the space-like trans- 
verse 4- vector e (e^ = —1) and hence decompose F into 
magnitude ao, envelope / = f{4') B^nd a constant tensor, 
f = n®e — e®n, 



F(<^) =ao/(0)f 



(28) 



with the dimensionlcss constant 4-vcctor n = kc/flo 
obeying — n ■ e — and n ■ u — n ■ uq ~ c. As 
a result, the square of f is simply = n®n with all 
higher powers vanishing due to = 0. In order for ([^5]) 
to be consistent with (fTS]) and ((27)) the average (. . .) must 
be defined in such a way that / is normalised, (/^) ~ 1. 
Defining a dimensionlcss gauge potential, A = eA/mc^, 
the field strength becomes 



F : 



n<gi 



A' -A',. 



(29) 



the prime henceforth denoting the derivative with respect 
to invariant phase (j). Comparison with (|28[) finally yields 

A' = aofe and = -{A''^). 

A more realistic laser model is provided by Gaussian 
beams, i.e. solutions of the wave equation in paraxial 
approximation [26j . The corresponding fields have non- 
trivial longitudinal and transverse envelopes resulting in 
the appearance of longitudinal field components. For this 
reason, the fields are no longer null and the charge dy- 
namics becomes more complicated. Charged particle ve- 
locities and trajectories in such fields have to be obtained 
numerically (27|. 



IV. SOLVING THE LL EQUATION 

Unlike the LAD the LL equation is a fairly standard 
equation of motion being second order in time deriva- 
tives. Hence, it requires two integrations and initial con- 
ditions for velocity and position. For the purposes of the 
present discussion it will be sufficient to perform only the 
first integration for which we need to provide the initial 
4-velocity, u{0) = uq. 



A. Neglecting Radiation Reaction 

To set the stage for a later comparison we first briefly 
recall the solution without RR, i.e. of the Lorentz force 
equation of motion (|13|) . We note that for any function 
/ = /(t) of proper time we may trade derivatives ac- 
cording to 



/ = /'0 = f'n = fno 



(30) 



where, in the last step, we have used the conservation 
law In terms of the dimensionlcss field strength F 

the Lorentz equation takes on the particularly compact 
form 



Fm 



(31) 
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This is, of course, integrated by a matrix exponen- 
tial which truncates at second order due to nilpotency, 
F^^ = 0. Employing the parameterisation (^5]) the solu- 
tion becomes 

ul = ua — cflo/ie + (ao/ie • Uq + ^allfcj n , (32) 

where the subscript "L" stands for "Lorentz" . The func- 
tion Ii (0) is the pulse shape integral 



Crucially, however, the equation ((37|) for the longitudi- 
nal velocity component completely decouples and, being 
first-order, can be solved by straightforward quadrature. 



1 



l + Rh 



(38) 



assuming the initial condition CIq ~ 1 and defining the 
shape integral 



h 



d(pf{(p) 



(33) 



Upon inspection of the solution ([5^ we note the following 
features. The velocity decomposes into transverse and 
longitudinal contributions given by the second and third 
terms on the right, respectively. If the initial velocity is 
longitudinal (like for a head-on collision) we have e • uq = 
0. In this case, the longitudinal velocity is quadratic in 
ao while the transverse component is always linear. 

One may rewrite the solution (|32|) in terms of the gauge 
potential A = oo/ie defined in p9|) which results in the 
neat expression 



uq — cA + (A ■ uq 



\A^c)n 



(34) 



From this expression one easily identifies the additional 
conserved quantity e • (u 4- cA) [28| corresponding to the 
transverse canonical momentum. As stated above, one 
has A • mq = for both fixed target and colliding modes. 
The quadratic contributions in p4|) arc actually positive 
as A is space-like, < 0. 



B. Including Radiation Reaction 

Upon including RR we have to solve the full LL equa- 
tion Q which we write in dimensionless notation as 



tlv! = Fu -I- r(f]F' - w)u , 



(35) 



withw = (u,F^tt)/c^, cf. 
frequency. 



11 = i7/ilo ^ k ■ u/k ■ uq 



26|) and ((27)) , and a normalised 



(36) 



Clearly, the LL equation (j35p is nonlinear in the unknown 
u. Remarkably, though, there is an analytic solution for 
a plane wave background, F = F((/)) [H, 113, HI- Let us 
briefly review its main steps using our compact notation. 

Recall that the Lorentz equation (|T3l) entails the con- 
servation law (|I5[) which is just = 1. In contradistinc- 
tion, a non-vanishing RR force, Frr = tqVF, implies 
that f2 is no longer conserved, but rather 



-rw 



(37) 



This is possibly the most significant new feature: In the 
presence of RR the electron will see a continuously chang- 
ing laser frequency during its passage through the pulse. 



(39) 



cf. ((33|) . It is worth noting that the RR parameter R 
from pB]) appears at this stage. As R ^ tq we should 
actually use the ultimate expression in (j38p in keeping 
with our philosophy of neglecting terms of order Tq . 

In any case we would like to point out that is a 
particularly nice signature for RR as it differs from unity 
only when a substantial amount of RR is present. In more 
physical terms, Q ^ I signals symmetry breaking in the 
following sense. Together with the longitudinal velocity 
M° — ~ O, cf. (|2ip . the longitudinal momentum, = 
pO _p3^ ceases to be conserved. As a result RR induces a 
breakingof translational invariance in the conjugate null 
directioi{3, = + x^. 

Let us continue with the LL equation and the remain- 
ing velocity components. The crucial technical trick (l5j 
is to introduce a new 4- velocity v via 



= n-'^u = (i + i?/2)u. 



(40) 



the longitudinal component of which is again conserved, 
k ■ V = ilo- Using (fST]) it is straightforward to see that 
the LL equation for v simplifies to 



= {n~'^W + rF' + rf^-^F^)?; . 



(41) 



As we know Ct^^ as a function of </> from ([38]) . the system 
(|4ip is indeed linear and easily solved via exponentiation. 
To arrive at the solution we use the parameterisation (j28p 
and discard all terms of order Tq (or r^). Noting that 
Uq — vq the solution for v may be then be written as a 
correction to the Lorentz solution (f5^ in the following 
way. 



UL + rAv + 0{r'^) 



(42) 



The RR term is explicitly given by 

Aw = -cao/2ie + (ao-^2ie • u^) + a^Iihic + a^hc) n , 
with the new shape integral. 



(43) 



(44) 



Recall the scalar product p ■ x = p^x /2 + p x'+/2 — ■ a;j 
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Comparing (|43|) and ((32l) one can identify precisely the 
same vector structure, the analogous longitudinal and 
transverse terms guaranteeing k-v = k-vo = flo = const. 
Obviously, in the limit of no RR (r — >■ 0) one has O — 1 
and Av such that ([32| is readily recovered from ([42]) 
and 

As a final comment we note that the dependence on 
proper time r is recovered by integrating ([55]) . which 
gives 



R I d^hiv) 



(45) 




Hence, in proper time r, RR leads to a phase shift com- 
pared to the Lorentz solution as r is no longer pro- 
portional to the phase, (j). 




FIG. 1. Laser pulse (|46|l for A'^ = 10 as a function of phase, 



AN ANALYTIC EXAMPLE 



With the analytic solution p3|) of the LL equation at 
hand we can readily analyse an example. The only re- 
maining technical difficulty is the evaluation of the pulse 
shape integrals ([33]), ([39]) and (|44p . It turns out that 
they may be performed analytically for the pulse shape 
function 



_ / sin4((/)/2Af)sin((/.) , 







< • 
else 



< 27riV 



(46) 



originally suggested in |30| (see also [3l|). The pulse 
(|46)) has a duration of (po = 2ttN (hence contains N cy- 
cles) and vanishes identically outside this interval. Thus, 
unlike a sine modulated Gaussian [29j, it has compact 
support (sec Fig.[T]). 



FIG. 2. Laser frequency Cl — k ■ u/k ■ uq as 'seen' by the 
electron during the passing of the pulse (A'^ = 10) as a function 
of phase, 0. Horizontal (black) line: Constant result psp 
without RR (no symmetry breaking). Decreasing (red) line: 
RR Solution (|38p signalling symmetry breaking. 



A. Symmetry breaking 

Let us first consider the behaviour of the symmetry 
breaking parameter Cl from ([55]). For this we need the 
integral I2, which is 




< < <^o 
else 



(47) 



The term 6I2 in (|47|) denotes a series of small amplitude 
sine functions which we do not explicitly display. All 
we need to know is that they vanish at the 'end' of the 
pulse, 6l2{4>o) = 0. Hence, inside the pulse = 1 — RI2 
drops linearly with small oscillations superimposed until 
it reaches a final plateau. For parameter values 70 = 100, 
ao = 150, ly = 10'^ (implying R = 0.022) and iV = 10 
the resulting behaviour is shown in Fig. ^ 

From (|47p the final plateau value, once the pulse has 
passed, is given by the simple expression 

35 

% = 1 - Rl2{M = 1 - S^riV— i?. (48) 

For the parameter values of Fig. [2] the numerical value 
for the plateau value is 51/ = 0.81. In general, assuming 
a head-on collision with 70 3> 1 one has VIq ~ 270, — 
2jf and the total relative energy loss becomes 



IA7I 
70 



1 - n 



f 



35 

256 



X 27rA^i? ~ NR 



(49) 



Apart from the numerical coefficient this is precisely our 
prediction ([^0)1 based on Larmor's formula. 

As an additional bonus ([48]) provides us with a crite- 
rion for when the LL approximation breaks down. This 
clearly is the case when NR becomes of order unity 
which, using ([TT]) . translates into 



10« 



(50) 
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FIG. 3. Laser frequency Q, as a, function of phase, (p, for a 
pulse with A'' = 10 and ao — 350. For tliis value, Cl becomes 
negative signalhng a breakdown of the LL approximation. 

For 7o :3> 1 one has e''" = 7o(l + Pq) ~ 270, so, when 
7o = 100 in Fig.[2]wc expect our approximations to break 
down when — 10^ or oq — 320. This is indeed borne 
out by Fig. H 

B. Varying initial conditions 

Let us finally look at the energy transfer dynamics in 
more detail. We want particularly to compare the two 
scenarios of fixed target and colliding modes which can 
both be described by the choice ([5]). All we have to do 
for a fixed target (electron initially at rest) is set /3o = 
and 7o = 1. We are interested in the electron energy as 
a function of 0, which describes its "history" during the 
passing of the pulse. This energy may be written as 

Eg = mcvP = ^mc? , (51) 

so we just have to monitor the behaviour of the electron 
gamma factor, 7(0) = vP{(j))/c. In the LL case this is 
governed by a nonlinear system, so we expect a significant 
dependence on initial conditions. This is indeed what 
happens. Let us first work out the analytical expressions. 

In analogy with (|42|) we split into a Lorentz and RR 
part, 

7 = 7L + A7 , (52) 

with the Lorentz contribution p2p always yielding an 
increase, 

7L=7o + ^e-«''a2/2>7o. (53) 
The leading RR correction is 

A7 = '^aval[h.f + alh + h{l - 7oe«") } , (54) 
with the new combination of shape integrals, 

h = h hl-\llh- (55) 
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FIG. 4. Electron gamma factor as a function of phase (j) in 
fixed target mode (70 = !)■ At the given vertical scale the 
curves with and without RR are indistinguishable. In both 
cases there is no net energy transfer. 

For the pulse ([46]) I4, turns out to be positive with com- 
pact support while Ii f can have either sign. All shape 
integrals are of order unity. Hence, for the case of in- 
terest (large ap ^ 1) the positive 1^ term dominates the 
Iif term. As a consequence, the sign of A7 is entirely 
determined by the last term, that is, by the initial condi- 
tions! The simplest case is the fixed target mode (FTM, 
7o = 1 = e*"") for which the crucial term vanishes and 
the RR correction is never negative, 

2 

A7FTM - 2a;^ao/4 > . (56) 

Asymptotically, once the pulse has passed {(f) > (po) we 
have A7FTM = as /4 has compact support. Hence, for 
FTM (and only for FTM!) there is no net energy transfer 
between electron and laser pulse (Fig. |4|). 

The situation is different in colliding mode (CM). As- 
suming 7o 3> 1 we have 

2 

A7CM ^ ^on^al (alh - 270/2) . (57) 

Unlike I4, I2 takes on a nonzero constant value after the 
pulse {4> > 0o). Hence, there is an asymptotic energy 
loss, 

A7CM/70 - -^Q:'^ao7o^2(0o) - -NR , (58) 

so that we recover ([^^ having identified its sign. The 
behaviour of 7 in CM is depicted in Fig. [5] both with and 
without RR. 

Let us finally trace the behaviour of the RR correction 
A7 during the crossover from FTM to CM, i.e. with in- 
creasing 7o. This is depicted in Fig. [H] The upper panel 
shows the FTM (70 = !)■ In agreement with the 
RR correction is always positive so that there is energy 
gain within the pulse. Compared to the amplitude 7 (see 
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FIG. 5. Electron gamma factor as a function of phase <j) in col- 
liding mode (70 = 100). Upper curve (green): Lorentz equa- 
tion (no RR). Lower curve (red): LL equation (RR) which 
implies net energy loss. 



Fig. (31) this is rather small, A7FTM ^ 7- This interme- 
diate energy gain increases with oq, cf. but no net 
energy transfer survives the passing of the pulse. 

If 7o is increased, 70 > 1, one enters CM (Fig. [SI central 
and lower panel). For intermediate values of 70 (central 
panel), A7CM stays positive as long as 4> is not too large. 
Once a sufficient number of cycles has passed, however, 
A7CM becomes negative which implies net energy loss. 

For sufficiently large 70, the RR correction A7CM is 
always negative, irrespective of the value of (f) (Fig. [6l 
lower panel). In summary. Fig. [6] thus shows the compe- 
tition of the second and third terms in (j54p . The latter is 
absent for 70 = 1 (FTM, upper panel). For intermediate 
7o (central panel) the ag term dominates for small 4> but, 
since I2 increases monotonically with 0, it overwhelms 
the sum for large cf) when /4 goes to zero. For sufficiently 
large 70 this holds for all (lower panel). 

This competition of ag and 70 has been noted before 
in the context of nonlinear Thomson scattering [32[ and 
also RR dynamics [s^l. There it was found that the 3- 
momcntum transfer changes sign at a "critical" value of 
ao — 270. For Thomson scattering this value defines 
an effective centre-of-mass system for which there is no 
energy transfer between electrons and laser photons. As 
a result, the theoretical emission spectrum degenerates 
into a pure line spectrum js^ . 



We close this section be reemphasising the crucial im- 
portance of working consistently to leading order in RR. 
Painful experience has told us that (wrongly) including 
higher order terms in the LL solution distorts the 
balance of terms in and one would not reach the 

conclusions above. 



A7 



Ay 








FIG. 6. RR correction to electron gamma factor as a function 
of phase <j) for ao = 150 but different values of the initial 
gamma factor. Upper panel: 70 = 1 (FTM). Central panel: 
70 = 40. Lower panel: 70 = 100. 



VI. CONCLUSION 

We have re-analysed the problem of radiation reaction 
by solving the Landau-Lifshitz equation analytically for 
an electron in an intense plane wave laser field. Such 
a field depends solely on the phase (j) = k ■ x 01, with 
the laser wave vector k pointing in z direction, on the 
light-front coordinate = ct — z. A particularly use- 
ful signature for radiation reaction is the laser frequency 



8 



as 'seen' by the electron, ft = k ■ u. This ceases to be 
conserved when radiation reaction is present and thus 
provides a clear signal for symmetry breaking: transla- 
tional invariance in the light- front coordinate = ct + z 
is lost. 

For a pulsed plane wave of finite duration in x~ (con- 
sisting of N cycles) the magnitude of the total change in 
51 (hence in longitudinal momentum, p^) obtained from 
the Landau-Lifshitz equation is well described by Lar- 
mor's formula for the radiated power. This has been 
corroborated by a careful study of the energy transfer 
dynamics and its dependence on initial conditions. The 
fixed target mode (electron at rest initially) is singled out 
as the only scenario for which there is no net energy trans- 
fer. Colliding modes with arbitrarily small initial velocity 
(70 > 1) always entail net energy loss. To arrive at these 
results it is of crucial importance to consistently trun- 



cate all expressions at leading order in radiation reaction. 
In this way it becomes obvious that the Landau-Lifshitz 
equation breaks down when jqUqN ~ 10^. Hence, if one 
wants to rely on the Landau-Lifshitz approximation, the 
electron energy 70, laser amplitude ao and pulse duration 
N cannot be increased arbitrarily. 

The integrability properties of both the Lorentz and 
Landau-Lifshitz equations seem quite intriguing. We 
plan to return to this topic elsewhere. 

ACKNOWLEDGMENTS 

The authors acknowledge stimulating discussions with 
S.S. Bulanov, A. di Piazza and A. Ildcrton. C.H. was 
supported by the Swedish Research Council Contract 
#2007-4422 and the European Research Council Con- 
tract #204059-QPQV 



H.A. Lorentz, La Theorie Electromagnetique de Maxwell [15 
et son Application aux Corps Mouvants, Arch. Neerl. 25, [16 
363-552 (1892), reprinted in Collected Papers (Martinus 
Nijhoff, The Hague, 1936), Vol. II, pp. 64-343. [17 
H.A. Lorentz, The Theory of Electrons, E.G. Teubner, 
Leipzig, 1906; reprinted by Dover Publications, New [18 
York, 1952 and Cosimo, New York, 2007. [19 
M. Abraham, Theorie der Elektrizitdt, Teubner, Leipzig, [20 
1905. 

P.A.M. Dirac, Proc. Roy. Soc. A 167, 148-169 (1938). [21 
H. Spohn, Dynamics of Charged Particles and their Ra- 
diation Field, Cambridge University Press, Cambridge, [22 
2004 [23 
F. Rohrlich, Classical Charged Particles, 3rd ed.. World 
Scientific, Singapore, 2007. [24 
K.T. McDonald, Limits on the Applicability of Classi- 
cal Electromagnetic Fields as Lnf erred from the Radiation 
Reaction, unpublished preprint, available from: [25 
http:/ /www. hep. princeton.edu/~mcdonald/examples/ [26 
L. D. Landau and E. M. Lifshitz, The Classical The- [27 
ory of Fields (Course of Theoretical Physics, Vol. 2), 
Butterworth-Heinemann, Oxford, 1987. [28 
J. G. Russo and P. K. Townsend, J. Phys. A 42, 445402 [29 
(2009) 

H. Spohn, The Critical manifold of the Lorentz-Dirac [30 
equation, Europhys. Lett. 49, 287 (2000) 

S. E. Gralla, A. I. Harte and R. M. Wald, A Rigorous [31 

Derivation of Electromagnetic Self-force, Phys. Rev. D 

80, 024031 (2009) [32 

C. Bamber, S. J. Boege, T. Koffas, T. Kotseroglou, 

A. C. Melissinos, D. D. Meyerhofer, D. A. Reis, W. Ragg [33 

et al, Phys. Rev. D60, 092004 (1999). 

T. Heinzl and A. Ilderton, Opt. Commun. 282, 1879 

(2009) 

J. Koga, T.Zh. Esirkepov and S.V. Bulanov, Phys. Plas- 
mas 12, 093106 (2005). 



A. Di Piazza, Lett. Math. Phys. 83, 305 (2008). 

A. Di Piazza, K. Z. Hatsagortsyan and C. H. Keitel, Phys. 

Rev. Lett. 102, 254802 (2009) 

Y. Hadad, L. Labun, J. Rafelski, N. Elkina, C. Klier, 
H. Ruhl, Phys. Rev. D82, 096012 (2010). 
F. Sauter, Z. Phys. 69, 742-764 (1931). 
V. Yanovsky et al.. Opt. Express 16, 2109 (2008). 
The Vulcan 10 Petawatt Project: 
http://www.clf.rl.ac.uk/New-l-Initiatives/14764.aspx 
The Extreme Light Infrastructure (ELI) project: 
http://www.extreme-light-infrastructure.eu 
T. Heinzl, Lect. Notes Phys. 572, 55-142 (2001). 
J. L. Synge, University of Toronto Applied Mathematics 
Series, No. 1 (Univ. of Toronto Press, 1935) 
H. Stephani, Relativity: An Introduction to Special and 
Ceneral Relativity, Cambridge University Press, Cam- 
bridge, 2004. 

W. Becker, Laser and Particle Beams 9, 603 (1991). 

L. W. Davis, Phys. Rev. A19, 1177 (1979). 

C. Harvey and M. Marklund, Radiation damping in 

pulsed Caussian beams, to appear 

H.M. Lai, Phys. Fluids 23, 2373 (1980). 

C. Harvey T. Heinzl, N. Iji, K. Langfeld, Phys. Rev. 

D83, 076013 (2011). 

F. Mackcnroth, A. Di Piazza, C. H. Keitel, Phys. Rev. 
Lett. 105, 063903 (2010). 

T. Heinzl, A. Ilderton, M. Marklund, Phys. Lett. B692, 
250 (2010). 

C. Harvey, T. Heinzl, A. Ilderton, Phys. Rev. A79, 
063407 (2009). 

A. Di Piazza, K. Z. Hatsagortsyan, C. H. Keitel, Phys. 
Rev. Lett. 102, 254802 (2009). 



